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ABSTRACT 

We present a systemic analysis of all of the stellar mass black hole binaries (confirmed & candidate) 
observed by the Swift observatory up to June 2010. The broad Swift bandpass enables a trace of disk 
evolution over an unprecedented range in flux and temperature. The final data sample consists of 
476 X-ray spectra containing greater than 100 counts, in the 0.6 - 10 keV band. This is the largest 
sample of high quality CCD spectra of accreting black holes published to date. In addition, strictly 
simultaneous data at optical/UV wavelengths are available for 255 (54%) of these observations. 

The data arc modelled with a combination of an accretion disk and a hard spectral component. 
For the hard component we consider both a simple power-law and a thermal Comptonization model. 
An accretion disk is detected at greater than the 5a confidence level in 61% of the observations. 
Lightcurves and color-color diagrams are constructed for each system. Hardness luminosity and disk 
fraction luminosity diagrams are constructed and are observed to be consistent with those typically 
observed by RXTE, noting the sensitivity below 2 keV provided by Swift. The observed spectra 
have an average luminosity of ^ 1% Eddington, though we are sensitive to accretion disks down to a 
luminosity of 10^"^ LEdd- Thus this is also the largest sample of such cool accretion disks studied to 
date. 

The accretion disk temperature displays two peaks consistent with the classical hard and soft spectral 
states, with a smaller number of disks distributed between these. The distribution of inner disk color 
radii is observed to be continuous regardless of which mode is used to fit the hard continua. Even 
before spectral hardening corrections are applied, there is no evidence for large truncation of the 
accretion disk in the hard state, with > 75% of the detected disks consistent with inner radii less than 
< 30 Rg, and with over 90% having radii less than ~ 100 Rg. Plots of the accretion disk inner radius 
versus hardness ratio reveal the disk radius to be decreasing as the spectrum hardens, i.e., enters the 
hard state. This is in contrast to expectations from the standard disk truncation paradigm and points 
towards a significant contribution from spectral hardening. Using the diskpn accretion disk model, 
we demonstrate that the observed apparent changes in disk inner radius are consistent with changes 
in the spectral hardening factor alone. 

The availability of simultaneous X-ray and optical/UV observations for a number of systems, allows 
us to examine the role of disk irradiation via a modified disk model with a variable emissivity profile 
(i.e., T(r) cx r~P). The data reveal irradiation of the accretion disk to be an important effect at all 
times. The accretion disk is found to dominate the UV emission irrespective of the assumed hard 
spectral component and the temperature profile reveals evidence for irradiation, even in softer states. 
Overall, we find the broadband soft state spectra to be consistent with an irradiated accretion disk 
plus a corona. 

Subject headings: accretion, accretion disks - black hole physics - X-rays: binaries 



1. INTRODUCTION 

Accretion onto a black hole is a fundamental astro- 
physical process capable of producing copious amounts 
of energy from radio - 7-ray frequencies. This radiation 
in turn will have a profound impact on its surroundings. 
In the local Universe this is illustrated by the discov- 
ery of large jet blown cavities around stellar mass black 
holes (Gallo ct al. 2004; Pakull ct al. 2010). When in- 
tegrated over the lifetime of the black hole, this energy 
will have a significant impact on its environment, e.g., 
feedback from supermassive black holes (SMBH) in the 
quiescent state (< 10~^ Lsdd) is required to reproduce 
the observed properties of galaxies (Croton ct al. 2006). 
The observation of tight correlations between the SMBH 
and properties of the host galaxy are further evidence 
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for the importance of feedback from the black hole, and 
understanding this will provide insight into the process 
of Galaxy formation and evolution (Ferrarese & Merritt 
2000; Gebhardt ct al. 2000; Gultckin ct al. 2009). At 
higher redshifts still, the accretion power of early HMXBs 
could have an impact on the formation of the earliest 
stars and galaxies (Mirabel ct al. 2011). 

There are over 500 accreting binaries currently known 
in the Milky Way, though this sample is dominated 
by systems containing neutron star accretors. Observa- 
tions of other galaxies have also revealed large popula- 
tions of accreting objects (Fabbiano 2006), which can be 
used to reveal information on the star formation history 
(Grimm ct al. 2003) and stellar content of these galax- 
ies (Gilfanov 2004). The X-ray luminosity function of 
the Milky Way can be used to constrain the population 
and evolution of X-ray emitting sources in the Galaxy 
(Grimm et al. 2002), which is found to be dominated by 
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accreting binary systems with compact object accretors 
(BH, NS & WD, Revnivtsev et al. 2011a,b). 

Detailed study of the accretion flow in Galactic X-ray 
binaries can provide constraints on the nature of ac- 
cretion in the presence of a large gravitational poten- 
tial, and the theory of general relativity itself (Psaltis 
2008). In particular, they provide a nearby window on 
black hole spin (Miller 2007; McChntock et al. 2011), 
the production mechanism of the ubiquitous relativis- 
tic jet emission (Markoff et al. 2001, 2005; Fender 200G), 
and even the nature of supernovae (SNe) and gamma- 
ray bursts (GRBs; Miller, Miller, & Reynolds 2011). To 
date, there are 20 confirmed stellar mass black holes 
(i.e., with dynamical constraints on the mass of the BH) 
and greater than 20 strong candidate systems, e.g., see 
McClintock & Remillard (2006). 

BH X-ray binaries (XRBs) are typically transient sys- 
tems that are observed to brighten from < 10~^ LEdd— ~ 
0.1 Lgdd on timescales of weeks to months before 
returning to a low luminosity state. The Galac- 
tic black hole binaries have been classified based on 
their behavior in the soft X-ray band (0.1 - 10 keV). 
This resulted in the definition of a number of appar- 
ently discrete spectral states (Tanaka & Shibazaki 1996; 
McClintock & Remillard 2006). When a black hole is 
accreting at an appreciable percentage of its Eddington 
luminosity (Lx > 10"'^ Lsdd, where Lsdd ~ 10^^ erg s^^ 
assuming a 10 Mq black hole), 2 primary spectral forms 
are observed. The first of these is the low-hard state, 
characterized by a power-law continuum in the 2-10 
keV band with a spectral index F < 1.8, and a spectral 
cut-off at energies ^ 100 keV. The RMS spectral variabil- 
ity is observed to be high in this state. The second state 
is the so-called high-soft state. The X-ray spectrum is 
dominated by emission in the soft band (< 2 keV) con- 
sistent with a thermal blackbody component, a softer 
low luminosity power-law component at higher energies 
(F > 2.2), and low RMS noise (Tanaka & Shibazaki 
1996; Grove et al. 1998). The launch of the RXTE mis- 
sion (Bradt, Rothschild, & Swank 1993) resulted in a re- 
finement in our understanding of these spectral states, 
and the transitions between them, which have resulted 
in the definition of a number of intermediate spectral 
states based on both the spectral state and the tim- 
ing properties of the source (Homan & Bclloni 2005; 
McClintock k Remillard 2006; van dcr Klis 2006). 

The broadband emission from accreting black holes 
is correlated with the behavior at X-ray wavelengths, 
as clearly illustrated by the observed phenomenology of 
the radio jet (Fender 2006). Studies of jets from ac- 
creting black holes on all mass scales have led to the 
realization that the accretion disk (inflow) and the jet 
(outflow) are closely related, i.e., 'fundamental plane 
of black hole activity' (Merloni et al. 2003; Falcke et al. 
2004) and 'disk-jet coupling' (Fender et al. 2004). In 
particular it has been noted that (i) the low-hard state 
(Lx ^ 0.05 LEdd ) is typically associated with the pres- 
ence of a compact quasi-steady jet (ii) radio /jet emission 
appears to cease when the system enters the high-soft 
state (Lx > 0.1 Lsdd)- This suggests that the accre- 
tion disk inflow is intimately involved in the process of 
launching a jet, e.g.. Fender (2006). 

At higher energies (> 10 keV), observations by Integral 
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Fig. 1. — Exposure time histogram for all of the systems in our 
sample, where the systems have been grouped in 0.5 ks bins. The 
exposure peaks at approximately 1 ks, with ~ 70% of the observa- 
tions having an exposure time less than 2 ks and 90% less than 4 
ks. The total exposure time is 940 ks. 



through the bulge monitoring program (Kuulkers et al. 
2007) and pointed observations are helping to elucidate 
the nature of the high energy X-ray component. Suzaku 
observations have demonstrated the importance of hav- 
ing broadband spectral coverage, where the broad band- 
pass (0.6 - 300 keV) has proven to be crucial to simulta- 
neously constrain the soft accretion disk component, the 
hard X-ray emission and the reflection features, i.e., Fe 
K line, Compton hump e.g., see Makishima et al. (2008); 
Reynolds et al. (2010a); Reynolds & Miller (2010b) for 
high energy spectra and Tomsick et al. (2009); Reis et al. 
(2011a,b) for constraints on the reflection component. 

1.1. Specific issues 

Understanding the relationship between the accretion 
disk inflow and the jet outflow is one of the outstand- 
ing goals of modern astrophysics. The anti-correlation 
between the accretion state of the black hole and the 
presence or absence of a relativistic jet indicates that 
understanding the geometry of the accretion flow in 
each of the primary accretion states will be key, i.e., 
the low-hard state and the high-soft state. The cur- 
rent paradigm for the conflguration of these accretion 
states is consistent with a hot geometrically thin opti- 
cally thick multi-temperature blackbody accretion disk 
(Shakiu'a ck" Siuiyaev 1973) extending to the innermost 
stable circular orbit (ISCO) in the soft state. As the 
accretion rate and hence disk temperature drop, the in- 
ner edge of the accretion disk recedes from the ISCO 
and the inner region is filled with a geometrically thick 
ADAF (advection dominated accretion flow, Esin c^t al. 
1997), which produces the observed hard spectral compo- 
nent that dominates the low-hard state X-ray spectrum 
(McClintock & Remillard 2006). 

However, a number of recent observations contradict 
this picture and reveal evidence for the persistent of 
the accretion disk at the ISCO in the low-hard state, 
e.g.. Miller, Homan, & Miniutti (2006a); MiUer et al. 
(2006b); RykoffetaL (2007); Tomsick et al. (2008); 
Wilkinson & Uttley (2009); Reis, Miller, & Fabian 
(2009); Reynolds et al. (2010a); Reynolds & Miller 
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Fig. 2. — Absorbed flux versus time for all of the black holes observed by Swift up to May 2010, see the symbols indicated on the right 
of the plot and Table 4. All fluxes are measured in the 0.6 - 10 keV band via the flux command in XSPEC. The outbursts of these black 
hole binaries are observed to exhibit a dynamic range of 10'' — 10^ over the 0.6 - 10 keV energy range. 



(2010b); Rcis, Fabian, & Miller (2010). These observa- 
tions suggest that the accretion disk is not the primary 
driver of the observed low-hard 4^ high-soft state tran- 
sition. Indeed, the latest observations suggest that the 
disk does not recede from the ISCO until luminosities of 
< 10-3 LEdd (Tomsick et al. 2009). A number of alter- 
native scenarios are consistent with the accretion disk 
remaining at the ISCO, e.g., the hard X-ray flux may 
originate in the relativistic jet (Markoff et al. 2005) or 
from a hot corona (Gierliuski et al. 1999; Merloni et al. 
2002), and it is the turning on/off of the jet or evolution 
of the corona that drives the observed state transitions. 

Questions also persist at longer wavelengths. Re- 
cent observations have highlighted a problem in our 
understanding of the UV emission mechanism in X- 
ray binaries. Standard theory explains the UV emis- 
sion as direct emission from the accretion disk and 
observations of the classical ly^^^ spectrum in early 
lUE observations supported this picture (C^licng et al. 
1992; Schradcr et al. 1993). However, it was soon 
realized that the large flux of X-rays present in an 
X-ray binary should modify the standard disk spec- 
trum, with the X-ray emission irradiating the disk, 
where it is subsequently reprocessed - disk reprocessing 
(e.g. van Paradijs & McClintock 1994; King & Ritter 
1998; Hynes et al. 1998; Dubus et al. 2001; Russell et al. 
2006). Over the last decade the picture has been 
further complicated by the realization that the ra- 
dio jet may contribute significantly across the entire 
spectral energy distribution, e.g., Markoff ct al. (2001); 
Russell ct al. (200G). Furthermore, simultaneous obser- 
vations at optical and X-ray wavelengths have revealed 
complex correlation functions, hinting at the presence of 
a non-thermal component in addition to the accretion 
disk emission, e.g., Kanbach et al. (2001); Durant ct al. 
(2008); Gandhi et al. (2008). 



Dunn ct al. (2010, 2011) has recently carried out a 
comprehensive study of black hole X-ray binaries us- 
ing over 13 years of archival RXTE observations. In 
Dunn et al. (2010), hardness intensity (HID) and disk 
fraction luminosity diagrams (DFLD) were constructed 
for 25 black holes. Dunn ct al. (2011) focused on the be- 
havior of the accretion disk component in each of the out- 
bursts studied in the previous paper. In particular, they 
analyzed the observed distribution of temperature versus 
luminosity for each outburst and found it to be consis- 
tent with LoiskfxTpjg]^, and with previous studies, which 
found the inner radius of the accretion disk to be constant 
in disk dominated states. They do note a deviation from 
this relation at the entry to/exit from the disk dominated 
state, which they suggest is consistent with a varying 
color correction factor (e.g., Shimura & Takahara 1995, 
see §4.3ii for further details). 

Despite the significant progress made over recent years, 
a number of unanswered questions remain. In this study, 
we focus our attention on the following issues: 

(i) Does the accretion disk recede from the ISCO dur- 
ing the soft to hard state transition? 

(ii) What is the nature of the accretion geometry in the 
low hard state, i.e.. Is the accretion disk truncated? 
What is the source of the observed hard X-ray flux, 
a jet or a corona? 

(iii) What produces the UV emission, does it originate 
from the accretion disk or is there a non-thermal 
contribution? 

The Swift observatory was designed to discover and pro- 
vide rapid multi-wavelength follow-up of gamma ray 
bursts (GRBs). These same attributes allow detailed 
study of many other time variable phenomena, in par- 
ticular accreting stellar mass black holes. Due to the 
rapid pointing capability and low pointing overhead Swift 
provides the first detailed CCD resolution monitoring 
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of black hole binaries in outburst. Combined with the 
strictly simultaneous observations at optical/UV wave- 
lengths provided by UVOT, Swift is the ideal platform 
to study the accretion process in detail. In compari- 
son to RXTE, the low energy X-ray cutoff allows us to 
study the accretion disk across the entire outburst cycle, 
whereas previously the cool disk would have been outside 
the RXTE low energy cutoff (- 3 keV). 

In this paper, we describe observations of black hole 
binaries undertaken with the Swift observatory prior to 
June 2010. In §2, we describe the observations and ex- 
traction of source spectra and lightcurves. We proceed 
to analyze the data in §3 and the results are presented 
in §4. In particular, we find that there is no evidence 
for significant truncation of the accretion disk, at least 
to luminosities > 10~^ LEdd- In fact, the observed small 
scale changes in the disk inner radius are found to be con- 
sistent with variations in the spectral hardening factor, 
which mimic a varying inner radius. In the broadband 
data, we find evidence for an irradiated accretion disk at 
all luminosities probed, but increasing as we enter the 
low-hard state. These results are discussed in the con- 
text of models for the accretion flow in §5, and finally 
our conclusions are presented in §6. 

2. OBSERVATIONS 

Our data sample consists of Swift observations of all 
known and suspected stellar mass black hole candidates 
observed through August 2010 (see Table 4). Data for 
each system was obtained via the HEASARC archive 
service. 

2.1. Swift 

The Swift observatory (Gehrcls ct al. 2004) is a highly 
versatile multi-wavelength platform, providing rapid 
pointing capabilities and simultaneous broad wavelength 
coverage from optical to hard X-ray energies. The Swift 
observatory contains three instruments: 

(i) Burst alert telescope (BAT, Barthelmy et al. 2005): 
A coded mask hard X-ray (> 10 keV) imaging telescope, 
covering the energy band from 15 - 150 keV. Low reso- 
lution spectra are also provided at higher fluxes. 

(ii) X-ray telescope (XRT, Burrows et al. 2005): A 3.5m 
focal length Wolter I type grazing incidence X-ray tele- 
scope, containing a 600 pix^ CCD (~ 2.35" pix~^, i.e., 
23 arcmin^), sensitive over the energy range 0.2 - 10 keV 
at the telescope focus, with an effective area of ^ 135, 
65, 20 cm~^ at 1.5, 6.4 and 8 keV respectively. Bright 
sources are typically problematic to observe with CCD 
based detectors (see §3.1.1). The XRT uses a number of 
different read-out modes to avoid/mitigate this problem, 
facilitating observations of sources with fluxes in excess of 
300 ct s~^. These may be divided into windowed timing 
(WT) and photon counting modes. In windowed timing 
mode only spectral information is provided as the central 
8 arc-minutes of the CCD is collapsed into a single spa- 
tial direction, providing 1.8 ms timing resolution. This is 
the primary mode used at high fluxes, with the observa- 
tions presented herein primarily utilizing WTW2 mode. 
Photon counting mode (PC) is the standard method of 
reading out a CCD providing both spectral and 2 dimen- 
sional spatial information but with a limited time reso- 
lution of only 2.5 s. The majority of photon counting 
mode observations herein are in PCW2 mode. 



TABLE 1 

Exclusion radius as a function of count rate 



Photon Counting Windowed Timing 



0.5 


- 1 ct s"^ 


2 pix 


150 - 200 ct 


3 pix 


1 - 


3 ct s-i 


5 pix 


200 - 300 ct s-l 


5 pix 


3 - 


6 ct 


7.5 pix 


300 - 500 ct 


10 pix 


6 - 


9 ct 


10 pix 


> 500 ct s-i 


15 pix 



Note. — Radius of the inner excluded region as function of 
count rate for all spectra considered herein, see also Romano et al. 
(2006); Rykoff et al. (2007). 



(Hi) Ultra-violet/ optical telescope (UVOT, Roming et al. 
2005): Simultaneous observations are obtained at op- 
tical/UV energies with the UVOT instrument (30 cm 
Ritchie Chretien telescope). The detector is a complex 
system consisting of a photodiode, micro-channel plates 
and a CCD (see the instrument paper for details). The 
CCD provides a 17 arcmin^ FoV and a frame time of 
10.8 ms. Coincidence loss occurs at high count rates, 
typically those in excess of 20 ct s^^, but this is not an 
issue for the optical/UV counterparts to the black hole 
binaries in our sample. Observations are normally ob- 
tained in only an single filter; however, there are some 
sources with multi-wavelength lightcurves. The systems 
in our sample contain data in one or more of the 6 UVOT 
photometric fihers V, B ,U, Wl, M2, W2 (A^ « 5000A, 
4500 A, 4000A, 2800 A, 2400A, 2200 A). 

3. ANALYSIS 

We aim to use the data sample described herein (see 
§3.5) to probe the physics of accretion onto black holes. 
To begin, we outline the steps taken to extract the final 
spectral and photometric data products. 

3.1. Reduction procedure: XRT 

All data reduction and analysis was performed us- 
ing the HEASOFT 6.9 suite, which includes ftools 6.9, 
SWIFT 3.5 and xspec 12.6.0k. The latest versions of 
the relevant Swift CALDB files are also used. 

Data were downloaded from the HEASARC archive 
using the coordinates listed in Table 4, where the search 
radius was restricted to 10'. All raw-data were repro- 
cessed using the xrtpipeline to ensure usage of the lat- 
est Swift CALDB when producing the relevant observation 
event files. Spectra were then extracted using xselect 
from circular and annular regions in the photon count- 
ing mode data, while rectangular regions were used to 
extract spectra from the windowed timing mode data. 
Exposure maps were generated for each observation using 
the xrtexpomap task. These are particularly important 
due to the presence of a number of dead columns near the 
center of the XRT CCD due to a likely micro-meteorite 
impact in May 2005 (Abbey et al. 2006). 

Background spectra were subsequently extracted from 
neighbouring regions. Ancillary response function files 
(.arf) were created with the xrtmkarf task, while the 
appropriate response matrix (.rmf) was sourced from the 
calibration database. These files were then grouped to- 
gether with grppha and exported to xspec for spectral 
fitting. 
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3.1.1. pile-up 

Due to the photon counting nature of modern X-ray 
CCDs, they are strongly susceptible to a phenomenon 
known as 'pile-up'. This occurs at high fluxes, where 
multiple photons can impact a detector pixel at approx- 
imately the same time (i.e., less than the detector frame 
time), and be recorded as a single higher energy event. 
This will produce an excess of hard photons in the mea- 
sured spectrum, which as a result will appear harder than 
is actually the case. Proper interpretation of the ob- 
served spectra requires a strategy to mitigate the adverse 
pile-up related effects. For further details see Miller et al. 
(2010) for a detailed discussion of pile- up in X-ray CCDs, 
including Chandra, XMM and Suzaku in addition to the 
Swift XRT. 

In our analysis, pile-up was corrected in the standard 
manner, i.e., by the extraction of annular regions cen- 
tered on the source position. The excluded region radius 
was chosen to be conservative and agrees with those used 
by other groups, e.g., Romano et al. (200G); Rykoff et al. 
(2007). The exclusion regions utilized are listed in Ta- 
ble 1 along with the corresponding count rates. In order 
to confirm the effectiveness of this strategy, we confirmed 
that the chosen exclusion was such that the best fit spec- 
tral parameters did not exhibit significant variance in 
comparison to the next largest exclusion region. The fi- 
nal spectra that we choose to analyze are those where 
the quantity of interest is robustly on the conservative 
side. 

3.2. Reduction procedure: UVOT 

The optical & UV photometry required minimal ad- 
ditional processing (Poole et al. 2008; Brcevekl et al. 
2010). We begin by utilizing the level II pipeline pro- 
cessed image file (sk.img). Images are aspect corrected 
via the uvotskycorr tool, where the images are regis- 
tered to the USNO-Bl catalogue (Monet et al. 2003). 
Source flux is extracted from an aperture of 5" radius and 
the background is extracted from a neighbouring source 
free sky position. The counterparts to the X-ray sources 
are easily identified, primarily due to the transient na- 
ture of the optical/UV counterpart in combination with 
the accurate X-ray position provided by the XRT. Source 
flux and magnitudes are calculated using uvotmaghist 
tool. These fluxes are converted to XSPEC readable '.pha' 
files with uvot2pha, which are then exported for spectral 
fitting. 

3.3. Spectral models 

The are a large number of models available in xspec 
for modelling the X-ray emission from black hole binary 
systems. As outlined in the introduction, the standard 
model for black hole X-ray binaries involves the combina- 
tion of soft emission (< 2 keV) from a thermal accretion 
disk and hard X-ray (> 2 keV) emission from power-law 
like component, typically attributed to Compton scat- 
tering of the thermal photons from the accretion disk to 
higher energies (McClintock Remillard 200G), or per- 
haps emission from a jet (Markoff et al. 2001, 200-5). In 
an effort to model the data in a manner that minimizes 
bias, we choose two complimentary approaches: 
(i) Phenomenological - Here we aim to determine the em- 
pirical relationships displayed in the data, which can be 
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Fig. 3. — Histogram of xi for our best fit pha(diskbb+po) model. 
The dasfied histogram is the model containing the power-law com- 
ponent alone, while the blue histogram denotes the best fit model 
after the addition of the diskbb component. An accretion disk is 
required at greater than the 5cr confidence level in > 6f% of the 
observations. Fits with alternative continuum components reveal 
similar results, see Table 2. 



used to inform/constrain theoretical models and charac- 
terise source types. For example, is the system spectrally 
hard or soft? A very soft X-ray spectrum is one of the 
characteristic signatures of an accreting black hole in the 
soft state. Additionally, we can create color-color and 
hardness intensity diagrams to search for characteristic 
variability and search for correlation at other frequencies, 
i.e., optical/UV. 

(ii) Physical - Using these models we will attempt 
to constrain the intrinsic properties of the accretion 
flow, e.g., accretion disk temperature and inner radius. 
For the soft component, we assume a standard steady 
state, geometrically thin, optically thick accretion disk 
(Shakura & Sunyaev 1973), described by the diskbb 
model (Mitsuda et al. 1984; Makishima, et al. 1986). Ex- 
tensions to this base model are also considered, e.g., cor- 
rections for GR effects in the inner disk region close to the 
black hole (diskpn, Gierlinski et al. 1999) and account- 
ing for possible changes in the temperature profile of the 
disk due to, for example, irradiation by the corona or 
flux generated otherwise close to the black hole (diskpbb, 
Mmeshige et al. 1994). 

The hard X-ray component, i.e., that which typically 
dominates above ~ 2 keV, is modelled using 2 distinct 
prescriptions. The most basic model for the hard X-ray 
flux is a simple power-law, characterized by a spectral in- 
dex, r. Alternatively, the hard component may be mod- 
elled assuming it originates via scattering of the low tem- 
perature seed photons from the accretion disk by a high 
temperature corona. Detailed Comptonization models 
including reflection compps (Poutancn & Svensson 1996) 
and non-thermal electron contributions eqpair (Coppi 
1999) are available; however, given the upper energy 
bound of the XRT (10 keV), we choose to use the most 
computationally efficient thermal Comptonization model 
(comptt, Titarchuk 1994). 

Where simultaneous UVOT/XRT data are available, 
joint fits are undertaken in an effort to constrain the na- 
ture of the broadband X-ray to opt/UV spectrum. The 
extinction in this case is modelled using the xspec phabs 
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Fig. 4. — Histogram of the observed luminosities for all of the 
observations in our sample, where the distances and masses listed 
in Table 4 are assumed. The unabsorbed luminosity is consistent 
with a normal distribution peaking at ~ 1% Eddington. 



& redden models for the X-ray and opt/UV data re- 
spectively. The standard Galactic extinction curve was 
assumed (CardcUi ct al. 1989) along with the gas to dust 
ratio from Prcdchl & Schmitt (1995), i.e., the reddening 
was set to E(B - V) = Nh/5.3 x lO^i cm'^. We note 
that in sightlines towards the Galactic center the column 
density is highly variable and as such the standard ex- 
tinction curve may not be appropriate for all sources in 
our sample. 

3.4. Analysis procedure 

For our baseline model, we utilize a blackbody ac- 
cretion disk in addition to a power-law to model the 
X-ray spectrum (i.e., diskbb+po). This model will 
also aid in the comparison of our results to the signif- 
icant archive of RXTE observations, e.g., (Dunn et al. 
2010). To model the absorption by intervening neu- 
tral hydrogen, we use phabs, where throughout this 
paper the abundances and cross-sections assumed are 
bcmc (Balucinska-Church & McCammon 1992) and angr 
(Anders & Grevesse 1989) respectively. The column den- 
sity is held fixed at a value consistent with previous ob- 
servations, see Table 4 for details. This is in agreement 
with the observed absence of significant local absorption 
in these systems (Milker ct al. 2009), with obvious excep- 
tions for the case of the wind-accreting HMXBs. 

Spectra with less than 100 total counts were ignored 
in all further spectral fitting. Such spectra are typically 
those obtained with the XRT operating in auto-mode, 
i.e., the first observation takes place in windowed tim- 
ing mode and if the count rate is safe the observation is 
stopped and switched to photon counting mode, some- 
times resulting in an under-exposed or piled-up spec- 
trum. The best fit statistic is chosen based on the to- 
tal number of counts in the spectrum. Standard chi- 
squared fitting is used for observations with more than 
500 counts. All spectra with sufficient counts to meet our 
'X^ fitting criteria were binned so as to have a minimum of 
20 counts in each bin via the grppha task thus ensuring 
that use of this statistic is valid. The Cash-statistic was 
utilized for those observations with less than 500 total 
counts (Cash et al. 1979). Finally, the spectral fitting is 



Fig. 5. — Histogram of the peak luminosity observed from each 
of the systems in our sample, where the distances and masses listed 
in Table 4 are assumed. Although the overall number is small, this 
distribution appears to peak above 10% Eddington with a signifi- 
cant tail to lower luminosities consistent with previous observations 
(Chen ct al. 1997; Dunn ct al. 2010). 



restricted to the energy range 0.6 - 10 keV. 

We then proceed to fit each spectrum as follows: (i) 
the spectrum is loaded in xspec, the number of counts 
are determined and if this is greater than 100 the fitting 
statistic is then chosen, (ii) the initial continuum model 
is fit to the data assuming a constant value for the in- 
terstellar column density as discussed previously, e.g., a 
power- law - pha*po. Best fit parameters and the associ- 
ated errors are determined in addition to the model flux 
(iii) an accretion disk component is then added to the 
model in the previous step, the best fit parameters of this 
new model are determined, e.g., pha* (diskbb+po) . (iv) 
an F-test is carried out to determine if the new 2 compo- 
nent model is a statistically significant improvement on 
the previous single component model, if the second com- 
ponent is required error scans are undertaken and fluxes 
are calculated. 

The above fitting procedure was then repeated, replac- 
ing the continuum components above with those indi- 
cated earlier, i.e., diskpn or diskpbb in place of the 
diskbb and/or comptt to replace the power-law compo- 
nent. A subset of these fits were also repeated with the 
column density as a free parameter to check the sensitiv- 
ity of our fitting methodology to this parameter. As the 
XRT bandpass is limited to energies below 10 keV, we are 
insensitive to the presence of a spectral cut-ofi^ at higher 
energies. As such, we are unable to simultaneously con- 
strain both the electron temperature (kTc) and optical 
depth (r) of the comptt model. We therefore choose to 
freeze the electron temperature to 50 keV in all fits using 
the comptt component unless otherwise indicated. This 
will result in the production of a power-law like spectrum 
in the 2-10 keV range. Additional fits were repeated at 
other electron temperatures to confirm the robustness of 
the fits to changes in this parameter. In the case of the 
opt/UV -f X-ray fits, the basic procedure was similar. 
The X-ray and UVOT data are loaded into xspec and 
the model is initialized at the relevant best fit model from 
the X-ray fits alone. The extinction at optical and ultra- 
violet wavelengths is modelled using the redden com- 
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TABLE 2 
Best fit model comparison 



Model 


Xl < 1-5 


xl < 2.0 


pha* (po) 


0.47 


0.58 


pha* (comptt) 


0.53 


0.68 


pha* (diskbb+comptt) 


0.83 


0.90 


pha* (diskbb+po) 


0.88 


0.93 


pha* (diskpn+po) 


0.85 


0.94 



Note. — Results of the best fit models for our sample, where a 
number of components are considered in addition to those plotted 
in Fig. 3. The final two columns indicate the percentage of systems 
with a best fit statistic less than indicated. 



ponent, i.e, pha(po)+redden(po) , and we relate this to 
the X-ray column density following PrcdoliI fc Schniitt 
(1995), i.e., E(B - V) = Nh/5.3 x lO^i cm-^. 

All of the above models were defined and fit using cus- 
tom tcl scripts. Errors where quoted are calculated via 
the error command and are equivalent to the 90% con- 
fidence interval unless otherwise explicitly stated. 

3.5. Sample characteristics 

Our sample comprises 27 stellar mass black holes. 
These range from those systems with meaningful con- 
straints on the mass of the compact object, Mx ^ 3 solar 
masses (e.g., Cyg X-1, GRO J1655-40), to a large num- 
ber of systems that are more correctly labeled as candi- 
date black holes. These are systems that have displayed 
spectral/photometric attributes consistent with those of 
the known black hole systems, e.g., soft spectrum, tim- 
ing properties, lack of pulsations or X-ray bursts, but 
lack a dynamical constraint on the mass of the compact 
object, e.g., 4U 1957-hll, Swift J1753.5-0127. The re- 
sults presented herein are necessarily biased by those ob- 
jects with the majority of observations, i.e., Cyg X-1, 
GX 339-4, XTE J1817-330, XTE J1752-223 and SWIFT 
J1753. 5-0127. However, these object are distributed rel- 
atively uniformly over the relevant binary parameters, 
e.g., HMXB vs LMXB, Porb etc., and hence do not overly 
bias our results. The known system parameters for each 
of the black hole binaries in our sample arc listed in Table 
4. 

A typical Swift pointing is short with an exposure time 
less than 2 ks. For the purposes of this study, we ignore 
all pointings for which the total number of counts de- 
tected by the XRT is less than 100. Thus our sample 
consists of the following: 

• 476 observations of 27 BHs with > 100 counts in 
the 0.6 - 10 keV band. 

• 255 observations for which we have simultaneous 
detections at optical/UV wavelengths, in at least 
one band (V, B, U, Wl, W2, M2). 

The lower number of optical/UV observations reflects 
the high column density through which some of these 
systems are observed, with many systems having Nh ^ 
10^^ cm~^. In Fig. 1, we plot the histogram of exposure 
times for all of the observation considered herein. The 
distribution peaks at an exposure time of ~ 1 ks, with > 
70% and 90% of the exposures being less than 2 ks and 4 
ks in length respectively. In Fig. 2, we plot the absorbed 
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Fig. 6. — Detailed zoom in on the outburst decay lightcurve for 
the best sampled systems, where to has been set to the time of 
maximum flux indicated by the shaded region (see also Fig. 2). 
Each system has been offset by an arbitrary amount for clarity. 



flux VS time for of all of the systems in our sample. We 
are sensitive to fluxes over a dynamic range of 10^. 

4. RESULTS 

The primary results from this study relate to the anal- 
ysis of the X-ray and UV spectra of each of the black hole 
systems in our sample. We briefly discuss the results of 
the X-ray spectral fits here. 

The spectra were fit with 2 families of models, one with 
a power-law representing the hard spectral component 
and the other with the power-law component replaced 
with a Comptonization model. In Fig. 3, we plot the xt 
distribution for the spectra when fit with a power-law 
only, and after the addition of a disk component, i.e., po 
— > diskbb+po. It is immediately apparent that a disk 
component is statistically required in a large number of 
observations. In total, we had 476 observations with > 
100 counts that were fit with the above model, and of 
these 445/476 (-- 94%) returned a best fit xl < 2. An 
accretion disk is required at greater than the 5(T confi- 
dence level, as measured by an ftest, in 282/445 (63%) 
of the observations with a best fit < 2. A further 
30 observations also required a disk, but in this case the 
"best" fit is stiU poor, i.e., xl ^ 2. 

In Tabic 2, we list the percentage of spectra below 
best-fit cut-offs of < 1-5 and 2 for our sample of spec- 
tral models. We find that fits with the hard component 
represented by the comptt spectral model return best 
fit values consistent with that returned by those mod- 
elled with a po component. This is to be expected, as 
we are unable to place constraints on the electron tem- 
perature of the Comptonizing corona due to the lack of 
data above 10 keV. Nonetheless, the accretion disk com- 
ponent is still strongly required in this case, though the 
disk properties are slightly different to those found with 
the power-law component (see §4.3). Likewise, replacing 
the diskbb component above with a diskpn model does 
not provide a statistically significant improvement in the 
best fit results. 

Before proceeding to a detailed discussion of the results 
of our spectral analysis, we first present the photometric 
properties of our sample. 
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(a) diskbb+po (b) Model tracks (c) diskbb+comptt 



Fig. 7. — Color-color plots for the systems in our sample, (a) Tracks traced assuming the best fit diskbb+po model, (b) The standard 
tracks traced by various model components, i.e., accretion disk, power-law and Comptonization. (c) Tracks traced assuming the best fit 
diskbb+comptt model. Fluxes have been measured from the best fit models with the cflux command in XSPEC. Only those observations 
that returned a best fit 2 are plotted. 



4.1. Light-curve morphology 

The morphology of a black hole binary X-ray lightcurve 
contains information which may be used to constrain the 
outburst, e.g., outburst energy, M etc. In Fig. 2, we plot 
the absorbed flux for each of the 476 pointed observations 
in our sample containing greater than 100 total counts. 
We are sensitive to outbursts over a dynamic range of 
greater than 4 orders of magnitude, and including the 
extra-galactic systems (NGC 300 X-1, ICIO X-1, M33 X- 
7), we are sensitive to accreting stellar mass black holes 
over 6 orders of magnitude in flux. 

In Fig. 4 the luminosity histogram of all observations 
in our sample is plotted, in both absorbed and unab- 
sorbed form. There is a clear peak in the distribution 
at a luminosity of ~ 1% Eddington^. Our limiting lu- 
minosity is ^ 10~^ IjEddi though there arc a number 
of isolated observations an order of magnitude fainter 
(due to GRO J1655-40, which experienced a large ampli- 
tude outburst, see Fig. 2). The unabsorbed luminosity 
equals or marginally exceeds the Eddington limit for a 
small sample of observations, which correspond to the 
known Eddington limited persistent sources, including 
GRS 1915-^105 (Vicrdayanti ct al. 2010) and the extra- 
galactic black hole IC 10 X-1, while the 2005 outburst 



TABLE 3 

LIGHTCURVE MORPHOLOGY 



System 


tb 


Tl 


T2 


Comment 


GRO J1655-40 


115 


20 


6 


Largest amplitude 


Swift J1842. 5-1124 




30 




Slow decay 


XTE J1817-330 


119 


32 


5 


Brightens ~ 60 days 










Const after ~ 130 days 


XTE J1752-223 




20 




Highly variable 


H1743-322 






4 


Fast decay 


Swift J1539. 2-6277 


26 


35 


7.5 


Slow decay for ^ 20 days 


GX 339-4 


15 


35 


7.5 


Flattens at ~ 30 days 


Cyg X-1 


5 


5 


120 


Const after ~ 5 days 


Swift J1753.5-0127 


28 


20 


40 


Const after ~ 80 days 


GRS 1915+105 






1 


Very fast decay 


Note. — Parameters of the decay lightcurvos displayed in Fig 



6. All times are measured in days elapsed from the time of maxi- 
mum fiux to. 



2 LEdd = 1-3 X 103S (Mx/10 Mq) erg s" 



from the transient black hole GRO J1655-40 also ap- 
proached this limit. In terms of flux, the distribution 
is broader owing to the varying distances of the detected 
sources. This peak corresponds to an unabsorbed flux 
of ~ 5 X 10~^ erg s~^ cm~^. The distribution of the 
maximum luminosity observed during each outburst is 
displayed in Fig. 5. A peak is observed at a luminos- 
ity of ~ 20% LEdd, this is consistent with the earlier 
observations (~ 0.2 Lgdd Chen ct al. 1997, ~ 0.12 Lsdd 
Dtmn ct al. 2010); however, small number statistics pre- 
clude any deeper analysis. 

In Fig. 6, we plot the lightcurves for 11 outbursts 
in our sample, where we have set to to be the time of 
observed outburst maximum. Due to the variable and 
transient nature of the black hole binaries in our sample, 
we note that our estimate of to is inherently uncertain for 
those sources without observations prior to the peak, i.e., 
XTE J1817-330, H1743-322, GRS 1915-H105. In addition 
to uncertainty in assignment of to , our interpretation of 
the lightcurves is also biased by the sampling and the 
duration of the outburst for which the system continued 
to be observed. This is particularly noteworthy for the 
cases of Swift J1842. 5-1124, Swift J1539.2-6277, H1743- 
322, GRS1915-t-105, and to a lesser extent GRS J1655- 
40. The lightcurves are characterized by fitting with a 
model assuming an exponential decay with timescale ri 
(Chen ct al. 1997). A break is added at late times (it, T2) 
where necessary, i.e.. 

The results of these fits are listed in Table 3. Our sam- 
ple of lightcurves includes 2 persistent systems that dis- 
play plateau type lightcurves (Cyg X-1, Swift J1753.5- 
0127). Three systems display an exponential decay type 
lightcurve followed by a break a later times (GRO J1655- 
40, XTE J1817-330, XTE J1752-223). Swift J1842.5- 
1124 & Swift J1539. 2-6277 also display exponential decay 
type lightcurves but with a characteristic decay timescale 
much shorter than the 3 previous systems. H1743-322 
was likely observed late in the decay phase of the out- 
burst as the lightcurve is consistent with that of GRO 
J1655-40 after the late time break. Finally, we display 
2 'outbursts' from GRS1915-I-105. These events display 
a much shorter timescale than that observed from the 
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Hardness Ratio [ 



2„kBV] 



Fig. 8. — Hardness luminosity diagram for the black holes in 
our sample. The hardness ratio is calculated from the ratio of the 
unabsorbed hard to soft X-ray fluxes, calculated from the best fit 
model pha(diskbb+po) . While there is a hint of the characteristic 
'q'-pattern visible in RXTE observations, e.g., Dunn et al. (2010), 
there is considerable scatter. 
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Powerlaw Fraction [ 0.6 - 1 keV ] 

Fig. 9. — Disk fraction luminosity diagram for the black holes in 
our sample. The dominance of hard state sources in our sample 
is clearly visible. Those systems with a power-law fraction less 
than 0.7 typically require the presence of a low temperature disk 
component. The spectra for which no disk component is required 
can be seen with a power-law fraction of 1. 



other systems, although the event is roughly consistent 
with that seen from Cyg X-1 during its transition to the 
soft state. Perhaps such flaring type variability is a char- 
acteristic of persistent accretion flows. 

The characteristic exponential decay timescales is sim- 
ilar across a number of systems, lying in the 20 - 40 day 
range. This is consistent with findings of the study of 
XRB lightcurves by Chen et al. (f 997). The typical out- 
burst durations of ~ 50 - 100 days (Fig. 6), imply a 
total outburst fluence in the range lO''^^ — 10'''^ erg s~^ 
for the galactic stellar mass black hole binary systems. 
Assuming the standard accretion efficiency (?y = 0.1), 
we find the typical amount of matter accreted per out- 
burst to be - 10-1° - 10"" Mq. Though the actual 
amount may be lower than this as the accretion effi- 
ciency is expected to be lower than 0.1 in the hard 
state, e.g., due to the presence of outflow/jets and/or 
an ADAF. Taking this into account, the measurement 
of an accretion rate onto the outer disk in quiescent 
black holes of Mx ~ 10-^° yr-^ would support this 
(McChntock et al. 1995; Orosz et al. 2011a). 

4.2. CC, HID, DFLD diagrams 

Color-color diagrams (Hasingcr & van der Klis 1989) 
are a useful means by which to characterise the spec- 
tral behavior and evolution of accreting black holes and 
neutron stars in a model independent manner. 

In Fig. 7, we plot color-color diagrams for the systems 
in our sample, based on the results of the diskbb+po 
model and the diskbb+comptt model. The soft colour is 
defined as (2.5 keV - 4.5 keV)/(0.6 keV - 2.5 keV), while 
the hard color is defined as (4.5 keV - 10.0 keV)/(2.5 keV 
- 4.5 keV). We plot unabsorbed fluxes as we want to 
avoid the detector dependent characteristics introduced 
by using the measured count rate. A number of our 
spectra contain some pile- up (§3.1.1), and correcting for 
this to obtain the true count rate is not straightforward. 
The colors are also impacted by the column density in the 
direction of each object, which varies considerably for the 
sources in out sample, i.e., Nh = 10^"'^ — 10^^ cm"^. The 
expected behavior for a number of different continuum 



components is plotted in the center panel, e.g., power- 
law, accretion disk etc. 

The Hardness luminosity diagram (HLD; Hoinan et al. 
2001; Belloni 2004; Fender et al. 2004) for the sources in 
our sample is plotted in Fig. 8. This diagram is lim- 
ited due to its reliance on count rates to calculate the 
hardness ratio. There will additionally be errors due to 
uncertainty in the column density, black hole mass, and 
distance estimates (Nh, Mx, d), which make it difficult 
to separate the different spectral states. If simultaneous 
timing information is available, then the HLD is very use- 
ful to characterise spectral evolution during an outburst. 
Nonetheless, we include it for comparison with previous 
work, e.g., Dunn et al. (2010). As we utilize a CCD de- 
tector, we choose not to calculate a hardness ratio from 
the spectral counts due to possible pile-up related effects. 
Hence, we define the hardness ration as the ratio of the 
hard to the soft flux, i.e., F2.0-10.0kev/F0.6-2.0keV- While 
the overall shape of the pattern traced by the objects 
is consistent with the 'q'-pattern observed with RXTE, 
there is considerable scatter, though this is to be ex- 
pected given the caveats outlined above and the differing 
bandpasses between the 2 instruments. 

In Fig. 9, we display the disk fraction luminosity dia- 
gram (DFLD) (Kalemci et al. 2004, 2006; Tonisick et al. 
2005; Kording et al. 2006) for the black holes in our Swift 
sample. The plot is generated from the best fit model 
diskbb+po over the spectral range 0.6 ~ 10.0 keV. Creat- 
ing this plot with the Comptonization model does not ap- 
preciably modify the resulting pattern traced by the ob- 
servations in the DFLD. The sample used by Dunn et al. 
(2010) is biased to higher luminosities in comparison to 
the current sample presented herein, which is dominated 
by hard state spectra. Dunn et al. (2010) calculated the 
luminosity in the 1 - 100 keV and 0.001 - 100 keV 
bands for the power-law and disk components respec- 
tively. This requires an extrapolation to energies well be- 
low the RXTE/PCA bandpass (i.e., below 3 keV), which 
will lead to errors as the accretion disk cools. Nonethe- 
less, the DFLD presented herein is qualitatively similar 
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(d) GX 339-4: diskbb+comptt 



Fig. 10. — Accretion disk properties for Cyg X-1 (left) & GX 339-4 (right) taken from the best fit diskbb+po (top) and diskbb+comptt 
(bottom) models. Displayed from top to bottom are the power-law index F (or optical depth, r), disk temperature Tin, change in the disk 
normalization relative to the minimum Adisijnomi) a-nd the hardness ratio defined as fo.6— 2.0kev/f2.0— lO.OkeV • 



to that measured by RXTE for a similar sample of black 
hole binaries Dium ot al. (201U). 

We can also ask if defining the hardness ratio in terms 
of energy ranges will mask changes in the source proper- 
ties. This is most likely an issue in the energy range 0.6 
- 3 keV, where there will be significant overlap between 
the hard (po, comptt) & soft (diskbb) spectral com- 
ponents. This has been investigated by comparing the 
number of soft /hard sources when wc define the hardness 
ratio as the ratio of the flux in each spectral component 
(fcomptt/fdisk) vs the the flux in the spectral bands above 
and below 2 keV. We find that for both definitions, the 
percentage of sources for which the above ratio is less 
than 1 is approximately 70%, i.e., both methods find ap- 
proximately 30% of the observations to be dominated by 
the soft component. For the hard sources, the maximum 
is observed at ^ 0.2 for those defined using fcomptt/fdisk- 

4.3. The accretion disk 

In this section, we focus on the 0.6 - 10 keV X-ray 
properties of the accretion disk as measured in our sam- 
ple, before also considering the available UV data. As 
outlined earlier (§3.5), an accretion disk is statistically 
required at greater than the 5a level in ~ 56% of our 
observations. We begin by considering disk models con- 
sistent with the standard Shakiu'a Sunyaev (1973) pre- 
scription. In Fig. 10, we plot the primary disk properties 



as a function of time for the best fit diskbb+po models in 
the case of the HMXB Cyg X-1 and the LMXB GX 339- 
4. These systems are chosen as examples of the accretion 
flow behaviors observed from a persistent and transient 
black hole binary system respectively. 

Cyg X-1 is a persistent X-ray source where the ^ 
15 M0 black hole accretes matter from the 0-type sec- 
ondary star (Orosz ct al. 2011b). It is normally detected 
at a flux of 0.2 ct cm'^ s-^ by Swift BAT in the 15 
- 50 keV band, while it is in the spectrally hard state. 
In Fig. 10, we see that Cyg X-1 was observed spar- 
ingly prior to the spectral softening which occurred in 
late 2010. The disk normalization (where the color ra- 
dius, rcoi oc norm*^'^) and spectral index are observed to 
increase as the spectrum softens while the temperature 
of the disk component remains fairly constant. For the 
Comptonization model the optical depth is observed to 
behave in the opposite sense to the power-law spectral 
index. We remind the reader that we have fixed the elec- 
tron temperature to be 50 keV due to the lack of data 
above 10 keV. Even though Cyg X-1 is a HMXB with a 
massive 0-type donor, in contrast to the the low mass 
main sequence stars typically detected in the transient 
black hole binaries, the basic X-ray continuum behav- 
ior described above is also observed from the persistent 
LMXB Swift J1753.5-0127. This suggests that the be- 
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(a) Disk temperature (b) Disk inner radius 

Fig. 11. — Left: Histogram of the accretion disk temperature for the disks detected in our sample. Right: Histogram of the accretion disk 
inner radius measured for the disks detected in our sample. The dashed vertical lines mark radii of 6, 10, 100 Rg assuming the distances 
listed in Table 4. The hard continuum component has been modelled in 2 separate ways (i) a power-law - po, and (ii) Comptonization - 
comptt. While both models are similar in the hard spectral region (> 2 keV), below this the Comptonization has a low energy roll over 
whereas the power-law does not. Nonetheless, there is considerable agreement between the derived disk properties in each case. 



havior we observe is driven by the inner accretion flow 
and not by the nature of the mass transfer, i.e., wind vs 
Roche lobe overflow. 

In contrast to the persistent systems, the transients 
as represented by GX 339-4 display discrete changes 
in measured disk temperature, consistent with the well 
known transitions between the low temperature hard 
state and the high temperature soft state, after bright- 
ening from the quiescent state (McClintock & Remillard 
2006). Since the launch of the Swift mission there have 
been 4 outbursts detected from GX 339-4 (Fig. 10 ~ 
right). Two of these outbursts are consistent with the 
system remaining in the hard-state for the duration of 
the outburst. Such hard-state outbursts have been de- 
tected from numerous black hole systems, e.g., GRO 
J0422-h32, XTE J1118-h480, GS 1354-64 among others 
(Brocksopp ct al. 2004). The other two outbursts are 
observed to show a clear transition to the accretion disk 
dominated soft-state. Here, we see that the normaliza- 
tion (color radius) of the accretion disk increases as the 
disk temperature decreases, though this effect essentially 
disappears when the hard component is modelled using 
the Comptonization model. The power-law index is ob- 
served to decrease in concert with the decreasing disk 
temperature, while again the optical depth behaves in 
the opposite sense, i.e., r increases as the disk tempera- 
ture decreases. 

The primary observational quantities providing us with 
information on the configuration of the accretion disk arc 
the inner disk color temperature (Tcoi) and the inner disk 
color radius (rcoi). The best fit values for the parame- 
ters are dependent on how the spectra arc modelled. In 
particular, the choice of model for the hard X-ray com- 
ponent can modify the measured disk flux, e.g., using 
a power-law (po) or thermal Comptonization (comptt ~ 
Titarchuk 1994) because the Comptonization cuts off at 
low energies consistent with the input seed photon tem- 
perature. We compare the 2 models above with the goal 
of determining, which model is a better description of 
the spectra. This will also allow us to probe the effect of 



incorrect modelling on the accretion paradigm (see §1.1). 
i. Disk temperature & radius 

In Fig. 11 (left panel), we plot the measured accre- 
tion disk temperature for both of the diskbb+po and 
diskbb+comptt continuum models. The primary dif- 
ferences are (i) that fewer observations require the ad- 
dition of an accretion disk component at greater than 
the 5cr significance level, i.e., 282 and 254 observations 
require a disk at the 5(7 level for the diskbb+po and 
diskbb+comptt respectively, and (ii) the disk component 
is typically cooler and has a slightly larger radius. 

In the diskbb+po fits there are a number of high 
temperature disks that are not required when the hard 
component is modelled using a Comptonization compo- 
nent. These high temperature disks arc seen from GRS 
1915+105, GX 339-4 and XTE J1752-223. In all of these 
cases, the spectrum is dominated by the hard compo- 
nent, i.e., these spectra appear to be consistent with the 
systems being in the very high state at the time of the 
observation (McChntock & RcmiUard 2006). 

The inner radius of the accretion disk may be deter- 
mined from the normalization of the diskbb model, i.e., 
norm = (rin[km]/d[10kpc])^cos6'. The disk radius thus 
measured is subject to a number of additional uncertain- 
ties, which bias the measured radius. The most impor- 
tant are (i) to account for spectral hardening, a mul- 
tiplicative correction factor of 1.7 is typically assumed 
(Shimura & Takaliara 1995, though sec below), and (ii) 
to account for the fact that the disk temperature does not 
peak at the inner radius (Makishima et al. 2000), which 
combined give us the corrected inner radius 



Rin[km\ ~ ^= — fliofcpc (2) 

V COS& 

In Fig. 11 (right panel), we plot the inner disk radius 
calculated for each of the accretion disks in the sample, 
assuming the distances listed in Table 4. The distribu- 
tion peaks at small inner radii, Rin ^ 10 Rg, with a tail 
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(a) dlskbb+po (b) diskbb+comptt 

Fig. 12. — Luminosity of the accretion disk component versus the disk temperature. Dashed Hues indicate the relation Lx cx: T**, as would 
be expected from a disk of constant emitting area. At large temperatures/fluxes this relation is followed by all of the systems with a small 
amount of scatter. At low luminosities, this is not the case with the disk temperature appearing to remain constant with changing flux. 
We note, at this disk temperature we are dominated by observations of the persistent systems Cyg X-1 and SWIFT J1753. 5-0127. 



to larger radii'^. Of particular note, the largest radii mea- 
sured in our sample are ~ 100 Rg. As a fraction of our 
sample, we find that 47% (24%) of the accretion disks 
have radii less than 6 Rg, increasing to 57% (36%) at 10 
Rg, 79% (74%) at 33 Rg, and 93% (91%) at 100 Rg for 
the diskbb+po and diskbb+comptt, models respectively. 

In Fig. 12, the relationship between the disk luminos- 
ity and temperature is plotted for both of the contin- 
uum models. If the emitting surface area of the accre- 
tion disk is constant (i.e., constant inner radius), then 
wc expect the systems to follow a relation of the form 
Ldisk oc T^. At higher temperatures and luminosities, a 
relation of this form (as indicated by the dashed lines) is 
observed; however, at the lowest disk temperatures (in- 
dicated by the shaded region), this no longer holds and 
instead the disk temperature is observed to remain ap- 
proximately constant while the luminosity changes by or- 
ders of magnitude. Modelling the accretion disk emission 
with the diskpn model (Gicrlinski et al. 1999), returns 
similar results. This model and deviations from the sim- 
ple Ldisk oc T"* relation are discussed in detail below. 

Examination of the temperature distribution in Fig. 1 1 
clearly reveals 2 distinct regimes, i.e., those at ^ 0.2 keV 
and those at a temperature of > 0.5 keV. These are con- 
sistent with the primary active accretion states, the hard 
state and the soft state respectively. We note that there 
is a deficit of systems in the temperature range between 
these two. The so called intermediate states are expected 
to reside here, i.e., those states during which the system 
transitions from the hard state to the soft state and vice 
versa (McClintock & Rcmillard 2006). In contrast, the 
inner disk radius distribution does not display obvious 
evidence for 2 or more distinct spectral states, instead 
being consistent with a single continuous distribution of 
low inner radii with a tail to higher radii, Rjn < 100 Rg. 
Finally, we note that this change in the disk tempera- 
ture and likely also the geometry /structure of the accre- 
tion flow takes place on timescales of days (e.g., see Fig. 
10), consistent with evolution of the inner accretion disk 



(Rin < 100 Rg) on the viscous timescale at those radii. 

The behavior of the hard X-ray components and their 
relationship to the observed behavior of the accretion 
disk was also examined, i.e., the power-law characterized 
by spectral index, F, and the comptonizing component 
characterized by the optical depth, r, see Fig. 13. Fo- 
cusing on the results of the diskbb+po fit, we find from 
left to right (a) disk temperature: a clear dichotomy ex- 
ists where we see a population of low temperature disks 
consistent with the hard state, and then a second popula- 
tion of hotter disks. Unfortunately due to our relatively 
narrow X-ray bandpass, the power-law is not well con- 
strained in those spectra where the disk is dominant, i.e., 
Tdisk 1 kcV. (b) The spectral index is observed to in- 
crease (i.e., soften) as the radius of the disk increases un- 
til a spectral index of approximately 2 is reached at which 
point the radius tends to be constant. The decrease of the 
accretion disk radius as the spectrum hardens points to- 
wards a changing spectral hardening factor, see §4.3. (c) 
At the lowest luminosities Lx < 10~^ LEdd, the spectral 
index is observed to remain constant, F ^ 1.8. At lumi- 
nosities above this the spectral index increases, though 
the rate at which it increases appears to vary system to 
system. Additionally, the spectral index may soften at 
the lowest luminosities; however, the number of spectra 
at these luminosities is too few to make any definitive 
statement. 

As noted previously, the diskbb+comptt model ex- 
hibits similar behavior to the diskbb+po model with the 
caveat that the optical depth, r, trends in the opposite 
sense to the power-law spectral index. This effect is strik- 
ingly apparent in Fig. 13. In the lower panels, we see 
that (d) the optical depth is observed to decrease as the 
disk temperature increases, (e) the inner radius of the 
accretion disk increases as the optical depth decreases, 
and (f) the optical depth of the corona decreases as the 
0.6 - 10 keV luminosity increases, i.e., as the spectrum 
becomes dominated by the soft thermal emission from 
the accretion disk. 



3 Rg = 14,8 (Mx/10 Mq) km 



ii. Spectral Hardening 
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(a) Tin vs r 



(b) Rin vs r 



(c) LEdd vs r 






(d) Ti„ vs T 



(e) Rin vs r 



(f) LEdd vs r 



Fig. 13. — Behavior of the hard X-ray component as characterized by the power-law index F (top row), and the comptonizing corona 
optical depth r (bottom row). The relationship between the above quantities and inner disk temperature, inner disk radius and Eddington 
luminosity arc plotted from left to right respectively. We note that the electron temperature for the corona has been frozen at 50 koV. 



In the previous section, the accretion disk emission is 
modelled as a standard multi-color blackbody accretion 
disk using the diskbb model. However, the emission 
from the accretion disk is not a true blackbody and in- 
stead will be modified as it scatters through the disk 
atmosphere. This emission may nonetheless be approx- 
imated by a blackbody. The temperature returned by 
the diskbb model. Tin, actually represents the tempera- 
ture of this modified blackbody which approximates the 
emergent disk spectrum. This temperature is related to 
the true effective disk temperature through the relation 

Tool = fcTeff (3) 

where fc is the spectral hardening factor 
(Shimura & Takahara 1995). In the standard sce- 
nario fc is assumed to be approximately 1.7 based on 
the work of Shimura fc Takahara (1995). Merloni et al, 
(2000) revisited this problem and found 1.7 < fcoi ^ 3, 
with the spectral hardening factor increasing as the 
spectrum became less disk dominated. Further theoret- 
ical study of the behavior of fc in the disk dominated 
state predicts fcoi ^ 1.7 for T < 1 keV (Davis et al. 
2005) and fcoi 2.2 for disk temperatures greater than 
1 keV (Davis et al. 2000; Done et al. 2008). With the 
exception of the study by Merloni et al. (2000), each 
of the other investigations focused on the soft state 
behavior of fc- 

In Fig. 14 (left), we plot the relationship between the 
hardness ratio and the inner disk radius, as calculated 
from the best-fit diskbb+po model. For the softer spec- 
tra (HR > 1), the radius is observed to be approximately 
constant for each source, whereas for the harder spec- 
tra (HR < 1), the radius is observed to decrease. In 
Fig. 14 (right), we plot the histogram showing the inner 



disk radius measured for each spectrum normalised to the 
smallest radius recorded from that system. The relative 
radius change is observed to peak at < 20x and lOx for 
the diskbb+po and diskbb+comptt models respectively. 

The apparent changes in the disk radius implied by the 
previous figures suggests that spectral hardening may 
play a role. In order to investigate this, the previ- 
ous continuum fits were repeated but the accretion disk 
model diskbb is replaced with diskpn (Gierliuski et al. 
1999). In this model, the inner radius is a free param- 
eter (which we fix at Rin = 6 Rg) and the normaliza- 
tion explicitly includes the spectral hardening factor, i.e., 
norm = (Mx/D)^cos(0)/f^. The results of these spectral 
fits, where a constant inner radius is explicitly assumed, 
are found to be consistent those using the variable radius 
diskbb model, see Table 2. Restricting ourselves to the 
predicted values for the spectral hardening factor men- 
tioned earlier (1.4 < /c < 3.0), we see that a change in 
the disk normalization by a factor of 10 - 20 is achievable 
through variations in the spectral hardening alone. 

Assuming that the radius variations observed in Fig. 
14 are solely due to variations in the spectral harden- 
ing factor, the diskpn model fits may be used to cal- 
culate fc- To do this, we assume that the minimum 
measured accretion disk normalization corresponds to 
the canonical spectral hardening factor, i.e., fc = 1.7 
(Gierlinski fc Done 2004). This minimum normaliza- 
tion corresponds to an accretion disk dominated observa- 
tion, and as the disk temperature/luminosity decreases 
we expect fc to increase as the continuum becomes 
dominated by emission from the hard X-ray component 
(Merloni et al. 2000). As the normalization cx we 
scale the other observations relative to this, hence al- 
lowing us to estimate the spectral hardening factor for 



14 



Reynolds et al. 



the observations in our sample. In Fig 15, we plot the 
histogram of calculated spectral hardening factors. The 
vast majority of observations are consistent with the the- 
oretically expected bounds on fc (1.4 < fc < 3.0), with a 
small tail of observations requiring a spectral hardening 
factor outside this range. 

In cases where the default value at normmin of fc ~ 
1.7 is incorrect, for example, in systems such as Swift 
J1753. 5-0127 that remained in the hard state at all times, 
a larger value of the spectral hardening factor may be 
warranted. This will not affect the observed pattern, 
which remains approximately the same; however, the am- 
plitude will be reduced. Therefore our choice of fc = 1.7 
at norminin is conservative, and the derived maximum 
hardening factors should be considered as upper limits 
in this case. In Fig. 16, we illustrate the behavior of 
the spectral hardening factor as a function of the disk 
temperature (left) and 0.6 - 10 keV luminosity (mid- 
dle). Finally, in the right hand panel we plot the cor- 
rected disk luminosity versus disk effective temperature, 
i.e., Lx (X T^g. 

4.4. UV emission 

i. Photometry 

In Fig. 17, we plot the relationship between the total 
unabsorbed X-ray luminosity in the 0.6 - 10 keV band 
and that detected in the optical/UV bandpasses. The 
differing number of points in each plot reflects the dif- 
fering observing strategies undertaken for each system, 
i.e., one or multiple UVOT filters. The optical emission 
is observed to be consistent with viscous emission from 
the accretion disk (Lopj/uv L^^"" '^), while in the UV 
we find evidence for both viscous emission and irradiated 
emission from the accretion disk (Lopt/uv L^;^). In the 
case of GX 339-4 (light blue), we find an approximately 
flat relation between the UV and the X-ray fluxes, partic- 
ularly apparent for luminosities less than ~ 10'^'' erg s^^. 
This hints at the presence of an additional source of flux 
during these observations, which we note are primarily 
from when the X-ray spectrum is power-law dominated, 
i.e., in the low- hard state. Neither the optical or the 
UV flux is consistent with a signiflcant contribution from 
the jet, which would be expected to produce a relation 
Lopt/uv oc L^^ (Russell et al. 2006). 

If instead we consider the ratio of the flux at X-ray 
to optical/UV- wavelengths, we find typical f opt/ fx 
10-5 - lO"*' in the V-band. The UV fiux ratio is simi- 
lar, though we caution that the red leak in the UVWl 
filter (Brecveld et al. 2011) may prohibit us from con- 
straining the true UV flux. Nonetheless there is evidence 
for differing behaviors. For example, SWIFT J1753.5- 
0127 is found to have a similar f opt/ fx, i-e., with in or- 
der of magnitude, during its entire outburst as observed 
by Swift thus far. We note though that this system is 
thought to have remained in the low-hard state since 
its discovery. In contrast, GX 339-4 displays a differing 
behavior, whereby when the system is in the soft state 
f opt/ fx ^ 10-5 but in the hard state f opt/ fx ~ 10"^ in 
the UVWl band. 

a. Broadband Spectroscopy 

An important question when modelling the UV emis- 
sion from an XRB is to determine whether the emis- 



sion is direct viscous dissipation from the accretion 
disk, reprocessed hard X-ray emission, or if it is 
a separate component entirely. Previous work has 
demonstrated that irradiation of the disk is impor- 
tant during the outburst of a black hole binary, e.g., 
van Paradijs & McChntock (1994); Dubus et al. (2001); 
Rykoff et al. (2007). Russell et al. (2006), using a large 
sample of black holes binaries, showed that the opt/nIR 
vs hard X-ray relationship was consistent with either an 
irradiated disk or jet origin. 

In order to investigate the relations observed above in 
more detail, the UVOT data were converted into xspec 
readable files and both the optical/UV and the X-ray 
data were fit simultaneously. To investigate the ori- 
gin of the UV emission, three models were considered 
(i) standard steady state disk plus a power-law compo- 
nent, i.e., diskbb+po, (ii) variable temperature accretion 
disk plus a power-law, i.e., diskpbb+po, and (iii) vari- 
able temperature accretion disk plus Comptonization, 
i.e., diskpbb+comptt. In all cases, we utilized phabs 
for the X-ray column & redden for the optical/UV col- 
umn assuming the standard Galactic extinction curve 
(Cardelli et al. 1989) and the dust to gas ratio from 
Prcdehl & Schmitt (1995). The column density was fixed 
at the standard X-ray value (Table 4). In those observa- 
tion where we had multiple opt/UV datapoints, it is pos- 
sible to constrain the column density. The typical best 
fit values were found to be consistent with the assumed 
X-ray column density within the errors in all cases. 

The results of the X-ray plus UV spectral fits may be 
summarised as follows: Model (i) diskbb+po: 267 broad- 
band spectra in total, and 196 of these are found to re- 
quire an accretion disk at greater than the 5(t confidence 
level. The number of broadband spectra which do not re- 
quire a disk component are 71, 60 and 46 at the 5cr, 4ct, 
and 3ct levels respectively. Model (ii) diskpbb+po: 264 
broadband spectra in total. Model (iii) diskpbb+comptt: 
257 broadband spectra in total. A similar number of 
models as in (i) above were found not to require an ac- 
cretion disk for model (ii). These spectra are discussed in 
detail below. A disk component is required in all of the 
diskpbb+comptt models as the corona cannot provide 
significant flux below the temperature of the input seed 
photons. When considering the quality of the model fits 
to the sample as a whole, they are found to be compara- 
ble, i.e., as measured by the respective xt distributions. 
However, the physical differences in the model compo- 
nents result in differing contributions from the accretion 
disk relative to the power-law, for example. 

As an illustrative test case, we plot the results of these 
fits to representative XRT-^UVOT spectra of GX 339-4. 
In Fig. 18; we plot fits to a soft state observation (obsid: 
00030919005). In fact, this observation is more accu- 
rately categorized as being a soft intermediate state, see 
McChntock & Remillard (2006). The unabsorbed model 
components are indicated by the thick grey dashed lines. 
All 3 models return a similar temperature for the accre- 
tion disk of Tin ^ 0.6 keV, and disk normalization of < 
1000, i.e, Rin < 40 dskpc O45 km. As can be seen in 
the figure, there is a clear difference between the model 
with the steady state accretion disk where the optical 
emission is dominated by the power-law component and 
those models with an irradiated disk, where the optical 
emission is dominated by the disk component. However, 



Swift black hole accretion survey 



15 



10000 




Hardness Ratio [ ° % o-,o ] jE,„ [ RJR„i„ ] 

Fig. 14. — Left: Normalization of the best fit accretion disk model (diskpn) versus hardness ratio, for the systems in our sample. The 
measured inner disk radii are low, consistent with Rjn < 100 Rg at all times. Right: Histogram of the ratio of the measured inner radius 
of the accretion disk in each observation relative to the smallest radius measured in any observation for each system. There are a small 
number of systems that appear to require a large change in disk radius; however, modelling the hard component with a Comptonization 
model instead of a simple power-law reveals these to be a systematic of the power-law model (see text). 



wc note here that the power-law index is poorly con- 
strained by the X-ray data, which is dominated by the 
disk emission. Both of the variable temperature profile 
models return non-standard profiles, i.e., T(r) cx r~*''^; 
however, the standard disk model provides a statistically 
comparable fit to the data. We must now consider the re- 
turned parameters and their implications. In particular, 
we note that both of the models that utilize a power-law 
to generate the hard X-ray emission, return hard pho- 
ton indices (F ^ 1.8). Such hard power-law emission is 
in conflict with the well known soft power-laws normally 
observed in the soft disk dominated states, i.e., F > 2.2 
(McClintock Rcmillard 2006). Hence for the soft state 
the data actually favour a model consisting of an irradi- 
ated accretion disk and a corona, i.e., diskpbb+comptt. 

In Fig. 19, we plot fits to a hard state observation of 
GX 339-4 (obsid: 00030943002), where the unabsorbed 
model components are again indicated by the thick grey 
dashed lines. In contrast to the soft state observation 
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Fig. 15. — Histogram of the color correction factor assuming 2 dif- 
fering hard X-ray components, i.e., a power-law and Comptoniza- 
tion. The dashed grey line indicates a color correction factor of 
fcol = 1.7. 



above, here the diskbb+po model returns a significantly 
poorer fit in comparison to the diskpbb models, i.e., 
xt ^ 1-3 vs 1.0. This is due to an inability of this 
model to re-produce the observed optical/UV slope as 
may be clearly seen in the residuals (middle panel) . The 
models with the diskpbb component are both observed 
to provide fits of a statistically similar quality. The re- 
turned accretion disk parameters are also similar, e.g.. 
Tin ~ 0.17 keV,Rin ^ 100 dgkpc ^'45 km. Crucially in 
both models, a non steady state accretion disk is re- 
quired. The measured temperature profile for the accre- 
tion disk is p = 0.52 ± 0.01. This temperature profile is 
consistent with that expected from an accretion disk that 
is being irradiated. Such irradiation is expected to occur 
in the hard state, when the X-ray emission is dominated 
by the hard X-ray component, i.e., E > 2 keV. Finally, 
we note the crucial role played by the broadband cover- 
age provided by the XRT plus UVOT combination. In 
particular, the coverage from optical to UV wavelengths 
is crucial in allowing us to constrain the temperature pro- 
file of the accretion disk, which in this case is consistent 
with an irradiated accretion disk, i.e., T cx r"" '^. 

The basic properties of the accretion disks (Tcoi, rcoi) 
in these fits are consistent with those measured in fits 
made to the X-ray data alone. However, it is important 
to note that it is impossible to constrain the temperature 
profile using the X-ray data alone, as the data will in- 
variably be consistent with a steady state disk within the 
statistical errors, e.g., Kubota et al. (2005). For the vari- 
able temperature profile model (T(r) cx r~P) the distribu- 
tion of the index differs between the Comptonization and 
power-law models reflecting the contribution by the hard 
component, in the power-law model case, to the opti- 
cal/UV emission. The average temperature profile (ilc) 
we measure is p ~ 0.59 ± 0.06 for the diskpbb+comptt 
model andp ^ 0.62±0.09 for the diskpbb+po model (see 
Fig. 20), though in both cases we arc biased by the large 
number of hard state observations in our sample as may 
be seen in Fig. 11. We emphasise, that irradiated disks 
are seen at all luminosities, with hard state observations 
typically requiring a more irradiated disk in comparison 
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Fig. 16. — Left: Color correction factor versus disk temperature. Middle: Color correction factor versus Eddington scaled luminosity 
for the systems in our sample. The Color correction has been calculated via the diskpn model assuming a constant inner disk radius of 6 
Rg. Right: Luminosity versus accretion disk temperature Toff, i.e., the colour corrected temperature of the accretion disk. The dashed 
lines denotes Lx oc T^. 



to soft state observations. 

There is no evidence for a dominant jet contribution 
at UV wavelengths as can be seen in Fig. 18 & 19. 
The power-law component is not dominant for any of 
the models in the hard state, and while it can dominate 
in the case of the diskbb+po model in the soft state, we 
do not consider this model to be valid as discussed above, 
i.e., the power-law would have to be much harder than is 
typically observed in a soft state. In all other cases the 
emission from the irradiated accretion disk dominates at 
UV wavelengths at all times. Though it is important to 
emphasise the fundamental difference implied by the two 
hard component choices. When the hard X-ray emission 
is most successfully modelled by a Comptonization com- 
ponent (i.e., soft state), then the UV emission is domi- 
nated by the accretion disc. In contrast, during harder 
states, where both the Comptonization and power-law 
components provide fits of comparable statistical qual- 
ity, two different physical scenarios arc suggested. For 
the Comptonization model, again all of the UV emission 
originates in the irradiated accretion disk, whereas for 
the power-law model the UV is dominated by emission 
from the irradiated accretion disk but with a contribu- 
tion from the power-law component that may be as large 
as 10% of the total continuum flux in the UV. Such a 
model may be consistent with timing observations that 
have found evidence for non-thermal flux at UV wave- 
lengths, e.g., Gandhi ct al. (2008). 

As stated earlier, there are a number of the broad- 
band spectra that arc consistent with a power-law alone, 
i.e, the addition of a disk component docs not signifl- 
cantly improve the flt. In the case of the diskpbb+po 
model 50/264 of the spectra are consistent with a power- 
law alone when we only mark an accretion disk as de- 
tected if it improves the fit statistic by greater than 
the 5(7 level as measured by an ftest. If instead we 
lower the criteria to a Scr improvement only 28 spectra 
are consistent with a power-law alone. The number of 
spectra consistent with a power-law alone is larger for 
the diskbb+po model; however, the excess relative to 
the diskpbb model is due to the inability of a standard 
steady state accretion disk to adequately reproduce the 
observed data. To examine these power-law spectra in 
more detail, we restrict ourselves to spectra with greater 
than 1000 counts and a best fit xt ^ 1-5, leaving us with 
24 spectra distributed across 6 systems (GRO J1655- 
40, GX 339-4, SWIFT J1753.5-1027, SWIFT J1842.5- 



1124, XTE J1752-223, XTE J1818-245). These systems 
broadly represent the contents of our sample, with no 
obvious biases towards system parameters, e.g., Porb- 
Though our knowledge of these parameters is limited in 
most of these systems, their observed behavior to date 
provides no evidence of any fundamental difference in 
comparison to the well known systems, i.e., GRO J1655- 
40. The average spectral index for this sample is found 
to be r = 1.66 ± 0.06 (Icr), while the average luminos- 
ity is (1.5 ± 0.7) X 10"'^ LEdd(lo'), both consistent with 
expectations for emission from a black hole in the low- 
hard state. The absence of a detectable accretion disk 
in these spectra may point toward a truncated accretion 
disk or conversely a disk that has cooled below our detec- 
tion limit. An alternative possibility, suggests that these 
spectra are dominated by emission form the broadband 
jet, e.g., Markoff ct al. (2001). 

5. DISCUSSION 

We have presented a comprehensive systematic analy- 
sis of all of the available Swift data on accreting stellar 
mass black hole binaries as of summer 2010. This is the 
most detailed study of accretion onto stellar mass black 
holes at luminosities of > 10"'^ LEdd- To summarize: 

• We analysed 476 CCD spectra of 27 systems, in the 
0.6-10 keV band. Approximately half of these had 
simultaneous coverage in the optical/UV band. 

• An accretion disk is required to successfully model 
the data in half of the X-ray observations. The ex- 
cellent sensitivity at low energies allows us to follow 
the evolution of the accretion disk to temperatures 
^0.2 kcV. This is the largest sample of such cool 
accretion disks collected to date. 

• The observed accretion disk components are con- 
sistent with maintaining inner radii < 100 Rg at all 
times, and typically less than 30 Rg, not accounting 
for any color corrections. More severe disk trunca- 
tion must take place at luminosities < 10^^ LEdd- 

• Independent spectral fits show that a constant in- 
ner disk radius but variable color correction factor 
is an equally good description of our sample. In this 
sense, the data are consistent with a disk remaining 
at (or close to) the ISCO for < lO^^ 

LEdd- 
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Fig. 17. — Left: Extinction corrected V-band versus X-ray flux. The optical data is consistent with expectations for viscous emission 
from the accretion disk (Lv— band c< Lx^""'"^). Right: Extinction corrected UVWl-band flux versus X-ray flux. The UV emission is 
displays evidence for irradiated disk emission (Ly-band <^ LSJ-^) at higher X-ray fluxes. 



• The simultaneous X-ray plus optical/UV data al- 
low us to constrain the temperature profile of the 
accretion disk. Irradiation of the accretion disk is 
important at all times, likely even at X-ray wave- 
lengths. 

• In the soft state, the broadband spectra are consis- 
tent with a spectrum dominated by emission from 
the accretion disk and an additional hard compo- 
nent consistent with Compton scattering of seed 
disk photons. 

• While an irradiated disk is preferred in the hard 
state, the data are unable to statistically dis- 
tinguish between power-law and Comptonization 
models for the hard flux. 

We discuss these results and their implications in detail 
below. 

5.1. Spectral hardening & disk truncation 

The accretion disk is the dominant spectral compo- 
nent at X-ray wavelengths, and the avenue through which 
matter is transferred from the donor star to the black 
hole. Improving our understanding of the characteris- 
tics and evolution of the accretion disk during an X-ray 
outburst is the primary motivation of this study. After 
all, the energy budget for any corona is set by the disk. 
The unique broadband data presented herein provide a 
new look at this problem. In particular, the low energy 
sensitivity of the XRT 0.5 keV) facilitates a detailed 
study of the accretion disk at low temperatures for the 
first time. This study contains the largest sample of such 
low temperature accretion disks observed to date. The 
addition of strictly simultaneous data at optical and UV 
wavelengths permits a detailed study of the temperature 
profllc of the disk and crucially its evolution as the disk 
temperature decreases. 

In total, an additional component consistent with a 
standard accretion disk is detected in over half of our 
observations. The excellent sensitivity of the Swift /XKT 
allows us to follow the evolution of these disks to tem- 
peratures ~ 0.2 keV (see Fig. 11). The disk evolution 



depends on how the hard X-ray emission is modelled, 
i.e., a simple power- law or Comptonization. This is most 
clearly illustrated in Fig. 10, where a prominent dis- 
crete change in the inner radius (rjn c>c -^norm) is ob- 
served when the hard X-ray emission is modelled with a 
power-law. In contrast, when a Comptonization model 
is used, the radius change is less pronounced. The tem- 
perature of the accretion disk follows a similar pattern. 
This behavior has been interpreted in terms of a model 
where, as the accretion rate decreases, the accretion disk 
recedes from the ISCO with the inner evacuated region 
fills with a corona/ADAF, and it is this re-configuration 
of the accretion flow geometry, which causes the classi- 
cal transitions between the hard and soft spectral states 
(Esin ct al. 1997). 

As discussed in §4.3, the accretion disk temperature 
measured from the observed X-ray spectrum is not the 
true effective temperature of the disk, but is instead the 
temperature of the disk modified as the emitted disk pho- 
tons pass through the disk atmosphere. This effect must 
be accounted for if we are to correctly interpret the data. 
This is done via the spectral hardening factor, fc, i.e., 
Tcoi = fcToff. Theory predicts the correction factor to 
be in the range 1.4 < fc < 3.0 (Shinuu'a k. Takahara 
1995; Merloni et al. 2000). 

In Fig. 11 (right panel), we show the inner disk radius 
calculated for each of the accretion disks in the sample. 
Regardless of the form of the hard spectral component, 
the inner disk radius is observed to be low at all times. 
The distribution of inner radii is observed to peak at 
R-in S 10 Rg, with a tail extending to larger radii. The 
largest radii measured in our sample are ~ 100 Rg, which 
rules out any large scale truncation of the accretion disk, 
for the luminosities probed in this survey (~ 10^'^ LEdd)- 
In Fig. 14, it is clear that the behavior of the apparent 
disk radius is consistent with spectral hardening: the ap- 
parent radius decreases as the source hardens. This runs 
counter to the prediction of ADAF models, which sug- 
gest that the disk should gradually recede in the hard 
state. Indeed, fits with the 'diskpn' model explicitly 
show that the spectra are consistent with a constant ra- 
dius while the correction factor changes in the manner 
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Fig. 18. — Simultaneous XRT & UVOT observation of GX 339-4 in the soft state, i.e., the accretion disk temperature is ~ 0.6 keV. The 
fits return of 1.23, 1.28 & 1.25 for models (a), (b)and (c) respectively. While the fits are statistically comparable, the power-law fit 
requires F ~ 1.8, in contrast to the expected value for a black hole in the soft state, i.e., F > 2. The models also attribute differing amounts 
of the optical/UV flux to the accretion disk component, ranging from 100% in model (a) to less than 10% in model (b). The preferred 
model (dlskpbb+comptt - a) requires an irradiated accretion disk, i.e., T{r) oc r"" ® (see text). The unabsorbed model components are 
indicated by the thick grey dashed lines. 




Fig. 19. — Simultaneous XRT & UVOT observation of GX 339-4 in the hard state, i.e., the accretion disk temperature is < 0.2 keV. The 
fits return xi of 1.02, 1.27,& 1.02 for models (a), (b) and (c) respectively. In this case the steady state accretion disk model in unable to 
reproduce the observed optical/UV spectral slope (panel b). The optical/UV emission is dominated by the accretion disk component in 
each model. The preferred models (those with the diskpbb component — a, c) both require an irradiated accretion disk, i.e., T(r) oc r~^'^ 
(see text). The unabsorbed model components are indicated by the thick grey dashed lines. 



predicted by Merloni et al. (2000). This finding does not 
rule out ADAF models in a broad sense; disk trunca- 
tion at much lower fractions of the Eddington limit is 
expected and may be partially confirmed by recent data, 
e.g., Tomsiclc ct al. (2009) on GX 339-4. 

For the first time we present a large sample of em- 
pirically determined color correction factors for a large 
sample of systems, see Fig. 15 & 16. The color correction 
factor is found to be consistent with the canonical value 
(fc ^ 1.7) for the majority of our observations, with a tail 
extending to larger values. When the hard X-ray compo- 
nent is modelled using the comptt model, the maximum 
value of the color correction factor is fc,max ^ 5. Figure 
16c plots the corrected disk temperature as a function of 
luminosity. At larger temperatures (larger disk fractions) 
the standard Lx c>c T^Jg is recovered; however, as we move 
toward lower disk temperatures (dominant hard compo- 
nent) a deviation is observed, hinting at an additional 
process here. I.e., is the disk finally beginning to recede 
as we approach these luminosities (Lx ^ 10"'^ LEdd e.g., 
see Tonisick et al. 2009). Future observations probing 
luminosities below 10~^ LEdd are necessary to provide 
further constraints on the evolution of the accretion flow 
geometry. 



The expected behavior of the colour correction factor 
in the low-hard state is uncertain. Merloni et al. (2000) 
showed that the color correction factor could increase rel- 
ative to the standard value of ~ 1.7. By how much is un- 
certain as in this regime the spectral continuum is neces- 
sarily dominated by emission from the hard component, 
i.e., corona or jet. In the standard soft state the X-ray 
spectrum is dominated by emission from the optically 
thick accretion disk, which is thought to extend down 
to the ISCO. In this spectral state the colour correction 
factor is expected to be roughly constant and observa- 
tions would appear to support this (Gierlinski & Done 
2004; Dunn et al. 2011). However, detailed theoretical 
treatments of emission from the disk in this state shows 
that evolution of the colour correction factor is likely 
(Davis et al. 2006; Done ct al. 2008). Gierhnski & Done 
(2004) analysed a sample of 10 black hole binaries and 
found the observed luminosity versus disk temperature 
relations to be consistent with a roughly constant spec- 
tral hardening factor for each system. The implied hard- 
ening factor was found to be 1.6 < fc ^ 2.3; although, 
we note that these observations were concentrated in the 
disk dominated state i.e., Tin ^ 0.5 keV. 

Spectral hardening appears to contribute to the ob- 
served disk properties at both high and low disk tern- 
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Fig. 20. — Histogram of the best fit accretion disk temperature 
profile for the dlskpbb model, i.e., T oc r~P. An irradiated tem- 
perature profile is required in a large number of observations. 



peraturcs/luminosities. At low luminosities, a dominant 
power-law component modifies the emitted disk spec- 
trum. At the highest luminosities, we are likely begin- 
ning to see the re-emergence of the hard component as 
the system moves to the very high state where in con- 
trast to the LHS the disk is now hot, i.e.. Tin ft ^■'^ keV 
(McClintock & Remillard 2006). In both cases, the color 
correction factor naturally increases as the non-thermal 
emission begins to dominate. At the lowest disk tem- 
peratures, we begin to approach the lower energy limit 
of the XRT, i.e., 0.6 keV (remember that the disk spec- 
trum peaks at ^ SkTjn). In the limit of low S/N some 
of the low temperature disks we detect could be system- 
atics caused by an excess of flux in some bins at the 
low energy end of our bandpass, which may mimic a low 
temperature disk component. Our decision to categorize 
a disk as detected only if it improves the fit by greater 
than 5(7 as measured by an ft est, should limit such er- 
rors in this sample of low temperature disks. Similarly, 
we have no reason to believe there is an issue with the 
XRT calibration; nonetheless, this effect requires further 
investigation, which we defer to a future study. 

Similar results have previously been presented by 
Dimn et al. (2011) based on an analysis of a large sample 
of RXTE spectra. The color correction factor was found 
to be approximately constant (~ 1.7) for disk dominated 
spectra, i.e., fdisk ^ 0.5. However, when studying the 
temperature luminosity relation for their sample, devia- 
tions from the expected behavior were observed at the 
beginning and end of disk dominated states, i.e., transi- 
tions from the soft state to the hard state and from the 
soft state to the very high state. Though the study was 
necessarily limited due to the low energy cut-off of RXTE 
> 3 keV, they were able to determine that changes in the 
color correction factor by at least Sfcoi 1 are possible. 

As noted above, there is no evidence that the accre- 
tion disk is truncated in the hard state, at least at the 
luminosities probed in this work (Lx ^ 10~^ LEdd)- 
This has long been known to be the case at the higher 
luminosities observed in the soft state, where the in- 
ner edge of the accretion disk is expected to lie at the 
ISCO, e.g.. Done Gierlinski (2003); Gicrlinski k. Done 
(2004). It is this constancy of radius that allows 



one to constrain the black hole spin via the accre- 
tion disk continuum emission, e.g., Steiner et al. (2010). 
In recent years evidence that the accretion disk is 
not truncated in the hard state has also been ac- 
cumulating, e.g., Miller, Homau, fc Miuiutti (2006a); 
Reis. Fabian. & Miller (2010). These results are in agree- 
ment with those presented herein, i.e., that the accretion 
disk in stellar black hole binaries is not truncated at least 
to luminosities greater than 10~^ LEdd- 

Previous efforts to model the behavior of the spectral 
hardening factor have primarily been focused on brighter 
more disk dominated states (Shimura & Takahara 1995; 
Merloni et al. 2000; Davis et al. 2005). The study by 
Merloni et al. (2000) is the only one that considers a sit- 
uation in which the disk does not dominate the X-ray 
emission and these simulations pointed towards an in- 
crease in the fc as the disk fraction decreased. Given the 
importance of this issue, more theoretical investigations, 
particularly in states where the hard component is dom- 
inant over the emission from the disk, are clearly war- 
ranted. The distribution of fc provided herein (see Fig. 
15 & 16) will provide a valuable guide to these studies. 

If the accretion disk does not truncate then we must re- 
consider the driving mechanism behind state transitions, 
in particular the hard and soft states, and how this is 
related to the production of a relativistic jet. The most 
promising means for powering these relativistic jets is the 
'Blandford-Znajek' mechanism, where the luminosity of 
the jet is Lbz oc B * r^ * f(a/M) (Blandford k Znajek 
1977). If the inner disk is present at the ISCO at all 
times as the above data implies, then this suggests that 
the accretion disk itself is not required to launch the 
jet. Indeed simulations have demonstrated the forma- 
tion of jets in the absence of a disk. For example, 
Barkov & Khangulyan (2011) analyze the formation of a 
jet launched from a black hole accreting from the stellar 
wind of a companion star. The driving mechanism be- 
hind the state transitions is less clear, though it is clear 
that a more detailed understanding of the evolution of, 
and interaction between, the hard component and the ac- 
cretion disk is key. The evidence presented herein points 
towards a non- varying disc geometry during the hard to 
soft state transition, as such, attention must turn to the 
corona and its coupling to the accretion disk as a function 
of mass accretion rate. 

5.2. Accretion disk irradiation 

The large non-thermal luminosities exhibited by ac- 
creting black holes leads one to naturally suspect this 
flux will modify the expected accretion disk spectrum. 
A standard steady state accretion disk will be observed 
to exhibit a radial temperature dependence of the form 
T(r) oc r-°-'^5 (Shakura & Sunyaev 1973). It is in- 
herently difficult to constrain the temperature profile 
of the accretion disk at X-ray wavelengths as the pri- 
mary modification to the accretion disk spectrum oc- 
curs at wavelengths to the red of the typical X-ray 
detector, i.e., UV/optical or E < 0.1 keV. This 
is illustrated in Kubota ct al. (2005) who studied a 
sample of ASCA observations of black holes in out- 
burst and found the temperature profile of the mea- 
sured spectra to be consistent with the standard T oc 
j,-o.75 accretion disk at all times, though there was evi- 
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dence, of low statistical significance, for a deviation to- 
wards a more irradiated disk in some of the observa- 
tions. In contrast, observations at optical and UV wave- 
lengths have revealed clear signatures for disk irradiation 
(Schrader ct al. 1994; van Paradijs <k McCIiutock 1994; 
RnsscU ct al. 2006; Rykoff ct al. 2007). Indeed, the ad- 
dition of X-ray irradiation of the accretion disk to the 
disk instability model is required in order to reproduce 
the observed outburst lightcurves of black holes XRBs 
(Dubus et al. 2001; Lasota 2001). 

In this work, we have utilized the combination of X- 
ray and UV/optical data provided by Swift to investi- 
gate the irradiation phenomenon in black holes binaries. 
The data may be divided between the 2 dominant accre- 
tion states (i) soft state: The X-ray luminosity is large 
(Lx ^ 0.1 LEdd) and the X-ray spectrum is dominated 
by the accretion disk, whose inner edge is consistent with 
lying at the ISCO. In Fig. 18, we plot fits to a soft state 
observation of GX 339-4. In contrast to the X-ray fit 
alone, where the temperature profile of the disk is con- 
sistent with the steady state expectation, the addition 
of the UV data reveals evidence for an mildly irradiated 
disk, i.e., T(r) oc r""-^. (ii) hard state: The X-ray lumi- 
nosity is typically lower (Lx ^ 0.1 L^dd) and the X-ray 
spectrum is dominated by a hard spectral component 
plus a cool accretion disk, where the precise location of 
the inner disk edge being uncertain. In Fig. 19, we plot 
fits to a hard state observation of GX 339-4. In this case 
the accretion disk is clearly irradiated, i.e., T(r) cx r""'^. 
The change in the temperature profile relative to the soft 
state is consistent with an increase in the irradiating hard 
X-ray flux, as expected in the hard state. 

In Fig. 20, we plot the distribution of accretion disk 
temperature profile obtained from the diskpbb broad- 
band spectral fits. There is a strong preference for a 
non-standard temperature profile irrespective of how the 
hard spectral component is modelled, i.e., po or comptt. 
There is no obvious evidence for a significant contribu- 
tion from a jet in the observations presented herein. In 
the hard state both the comptt/corona and po/jet pro- 
vide acceptable fits to the data. Fundamentally, we are 
unable to differentiate between these 2 physical scenar- 
ios with the available data, though appealing to previous 
timing observations (e.g., Gandhi et al. (2008)) leads us 
to favor the power-law model. In contrast, in the soft 
state the spectra do appear to rule out a power-law ori- 
gin for the observed hard X-ray component as this results 
in a much harder spectral index than is normally mea- 
sured in the soft state. This is in agreement with the well 
known behavior of relativistic jets, whereby they are ob- 
served to be present in the hard state, but appear to be 
absent in the soft state (Fcnidcr 2006). 

In light of these results, it may be necessary to re- 
consider the effect of the assumed temperature profile on 
the measured black hole spin, as even though the X-ray 
spectrum may be consistent with a standard steady state 
disk, when one considers the broadband spectrum devia- 
tions from this profile become apparent. This may mate- 
rially impact the derived disk parameters. For example, 
a single observation of the black hole binary 4U 1957-1-11 
is included in this study. When the broadband data is 
fit with the variable temperature profile disk model, the 
best fit temperature profile is consistent with an irradi- 
ated accretion disk with a lower disk temperature in com- 



parison to the steady state case. Recently, Nowak et al. 
(2011) have measured the spin of the black hole in this 
system via the continuum method using Suzaku spectra, 
where the spin was consistent with a near maximal value. 
Clearly, failure to correctly account for disk irradiation 
could bias the measured spin. We defer detailed consid- 
eration of this to a future study. 

Finally, we would like to re-emphasise the importance 
of the broadband coverage provided by Swift in con- 
straining the properties of the observed accretion disk 
emission, which is key if we are to understand the re- 
lationship between the accretion disk inflow and the jet 
outflow. 

6. CONCLUSIONS 

We have presented a comprehensive analysis of all of 
the data on stellar mass black binaries in the Swift public 
archive as of June 2010. This is the largest collection of 
CCD spectra of accreting stellar mass black holes pub- 
lished to date, amounting to 476 spectra of 27 binary 
systems. The excellent sensitivity and broad soft X-ray 
bandpass of Swift allows us to study the evolution of the 
accretion disk down to temperatures as low as ~ 0.2 keV 
for the first time. Approximately half of the spectra are 
found to require the addition of an accretion disk com- 
ponent at greater than the ba. 

There is no compelling evidence in the data presented 
herein, which samples the accretion disk evolution down 
to luminosities of ^ 10"'^ Lsddj for a large scale change in 
the inner disk radius during the transition from the hard 
to the soft state, i.e., Rin ^ 100 Rg at all times. Consid- 
eration of the evolution of the inner radius versus spectral 
hardness points towards spectral hardening being a fac- 
tor in the observed disk radius evolution. Taking this 
into account, we find that the observed apparent small 
scale inner radius changes can be accounted for when 
plausible spectral hardening corrections are applied to 
the data. 

Irradiation is found to be present in the majority of ob- 
served spectra. This irradiated disk component is found 
to dominate the optical/UV emission in both hard and 
soft states. The soft state spectrum is found to be con- 
sistent with a combination of an irradiated accretion disk 
plus corona, whereas for the hard state a significant con- 
tribution from a power-law component (i.e., a jet) cannot 
be ruled out at optical/UV energies. 

The data and analysis presented herein clearly demon- 
strates the value of the broadband observational capabil- 
ities of the Swift observatory, in the study of accretion 
physics. An analysis of the neutron star XRBs observed 
by Swift, and a comparison with the black hole binary 
systems presented herein will be presented in a compan- 
ion paper (Reynolds et al., 2012 in prep.). 
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bution made to the study of accretion physics by the 
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TABLE 4 

Confirmed black hole binaries: primary properties. 
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